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Abstract: The optical spatial solitons with ellipse-shaped spots have gen- 
erally been considered to be a result of either linear or nonlinear anisotropy. 
In this paper, we introduce a class of spiraling elliptic solitons in the 
nonlocal nonlinear media without both linear and nonlinear anisotropy. The 
spiraling elliptic solitons carry the orbital angular momentum, which plays 
a key role in the formation of such solitons, and are stable for any degree of 
nonlocality except the local case when the response function of the material 
is Gaussian function. The formation of such solitons can be attributable to 
the effective anisotropic diffraction (linear anisotropy) resulting from the 
orbital angular momentum. Our variational analytical result is confirmed by 
direct numerical simulation of the nonlocal nonlinear Schrodinger equation. 
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1. Introduction 

The nonlinear propagation of optical beams with ellipse-shaped spots has been discussed during 
recent years. The self-trapping beams with ellipse-shaped spots can be obtained by introducing 
either linear anisotropy or nonlinear anisotropy. Elliptic incoherent solitons have been reported 
in saturable nonlinear media 0~), in strongly nonlocal media |2 | and in photorefractive screening 
nonlinear media [3], where linear anisotropy comes from the anisotropic coherence function. 
Furthermore, it was predicted H that there exists an elliptical self-trapping beam for the ex- 
traordinary light in uniaxial crystals due to the anisotropic diffraction (linear anisotropy) (5) 16). 
On the other hand, coherent elliptic strongly nonlocal solitons were observed experimentally 
in lead glass Q where nonlinear anisotropy is achieved by rectangular boundaries in the trans- 
verse, and were also simultaneously and independently predicted when the nonlinear response 
function of the medium was assumed to be anisotropic [8|. And elliptical discrete solitons can 
form in an optically induced two-dimensional photonic lattice where the nonlinear anisotropy 
comes of enhanced photorefractive anisotropy and nonlocality under a nonconventional bias 
condition j9). 

Since optical solitons are the result of the exact balance between linearity and nonlinearity, 
the elliptic solitons can generally not exist in the media with both linear and nonlinear isotropy, 
and the ellipse-shaped beams always undergo significant oscillations in the propagation direc- 
tiom in such media iTlOl fTTI IT21 . It was predicted very recently 031, however, that the elliptic 
solitons with the initial orbital angular momentum (OAM) can exist in such media, and they 
will rotate along the propagate distance. Such solitons are unstable in the (local) cubic nonlin- 
ear media, but they can propagate stably in the saturable nonlinear media because the saturable 
nonlinearity can arrest the collapse instability lfl4l . 

Apart from the saturable nonlinearity, there is another mechanism that is nonlocal nonlinear- 
ity can arrest the catastrophic collapse of the self-trapping beams 1151 , So such class of spiraling 
elliptic solitons might also exis t in nonlocal nonlinear media, which will be confirmed theoret- 
ically in this paper. 

2. The variational solution of the spiraling elliptic soliton 

The propagation of optical beams in nonlocal cubic nonlinear media can be modeled by the 
following nonlocal nonlinear Schrodinger equation(NNLSE) lfl6l . 

2^ + + 0+2^Ayy^-^r,-T,O|A(r J nOI 2 dW = O, (1) 



where A(§,TJ,Q is a paraxial beam, R is the response function of the medium, £ is the lon- 
gitudinal coordinate, £, and 77 are the transverse coordinates, k = cohq/c is the wavenumber in 
the media without nonlinearity, yiq is the linear refractive index of the media, n% is the non- 
linear index coefficient. Through the dimensionless transformation x = £, /wo,y = r)/wa,z = 
£/(&Wq), y/ = Akwo(ri2/no) 1 / 2 ,R = WqR, where wo is the initial width of the optical beam, 
Eq.([T]) is expressed as in the dimensionless form 

.dw 1 / d 2 w d 2 w\ 

+ o l^T + ^T )+Anx l r=0, (2) 



dz 2 \ dx 2 dy 2 

where An = // R(r — r')| l/f(r') | 2 d 2 r' is the nonlinear perturbation of refraction index with 
r = xt x +yt y . R(r) — \j (2nw 2 n ) exp[— r 2 /2w 2 } is assumed in this paper, where w m is the char- 
acteristic length of the response function in the dimensionless system. 

The Lagrangian of Eq.© can be expressed as flTl L = 1/2 ff{y* dy/dz — 
y/dxi/* /dz)dxdy - H, where H is the Hamiltonian of this system, H = l/2 ff(\d\j//dx\ 2 + 
\dy//dy\ 2 — An I y/| 2 )dxdy. We introduce a trial function lfl3l . 
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where the Gaussian envelope is G{t ) = exp(— 1 2 /2), the phase is <j) = B(z)X 2 + &(z)XY + 
Q(z)Y 2 + (p(z),X =xcosp(z)+ysinp(z),Y = -xsmp(z) +.ycosj3(z) andPisthe power,P = 
// I \j/\ 2 dxdy. We can obtain the orbital angular momentum(OAM), M = lmJJ \jf* (r x V y/)d 2 r 
I /2P(b 2 — c 2 )®. Inserting the Gaussian ansatz ([3} into the Lagrangian, L can be analytically 
determined. Then using the variational approach, we can obtain that P' = 0, H' = 0, M' = 0, b' = 
2bB,c' = 2cQ,p' = (b 2 + c 2 )®/(b 2 -c 2 ) and 

m /_ b 2 + c 2 P[6b 2 c 2 + 5 W l(b 2 + c 2 )+4 W i} 

9 ~ 2b 2 c 2 + 87t[(b 2 + wl)(c 2 + wl)]y 2 ' 1 ^ 

where the primes indicate derivatives with respect to the variable z- So it can be found that the 
power, the Hamiltonian and the OAM of the system are conservative. We can determine the 
Hamiltonian of the system, H = P/A{b' 2 + c' 2 + n) , where 

n _ J_ ]_ 4b 2 a 2 Ac 2 a 2 P 

~ b 2 + C 2 + (ft2 _ c 2f + (pi _ c 2f ~ % ^r + w 2 n ) ( c 2 + w 2 m ) ' 

with o = M/P = 1 /2(b 2 - c 2 )©. 

Solitons can be found as the extrema of the potential Yl(b,c). Assuming b > c without loss 
of generality and letting dll/db = and dll/dc = 0, we can obtain the critical power and the 
critical OAM 

2^(l+p 2 ) 3 [(l+g 2 )(l+g 2 p 2 )] 3 / 2 (p 2 -l) 2 [l+g 2 (l+p 2 )] 1 / 2 
p[l + (6 + 45 2 )p 2 + (l+45 2 )p 4 ] ' e 2p[l + (6+45 2 )p 2 + (l+45 2 )p 4 ]V2' K > 

and M c = P C G C , where p = b/c and 8 = w m /b represent the ellipticity of the elliptic beam 
and the degree of nonlocality, respectively. The larger is 8, the stronger is the degree of non- 
locality. When P = P C ,(T = a c , the optical beam can propagate keeping its elliptic profile 
changeless and rotating stably. We can also obtain the rotation velocity D. c = j3' = 2(b 2 + 
c 2 )cr c /(fc 2 -c 2 ) 2 . When the semi-axes b and c are given, the critical power and the critical 



OAM of the spiraling elliptic solitons can be determined by Eq.©. One example is shown 
in FigOIa) with P c = 1.27 x 10 5 ,ct c = 0.561, C = 2a c /(b 2 - c 2 ) = 0.374 and Q. c = 0.623 
when b = 2.0, c = 1.0, w m = 15.0. Comparing two half widths obtained from variational so- 
lution, w x = (fe 2 cos 2 n c z + c 2 sin 2 £2 c z) 1 / 2 and w y = (c 2 cos 2 £2 c z + ^ 2 sin 2 i2 6 z) 1 / 2 , with those 
from the numerical simulation of Eq.© by using y(x,y,0) = [P c /(7tbc)] l / 2 exp[—x 2 /(2b 2 ) — 
y 2 1 (2c 2 )] exp(/0 c xy) as the input beam at z = 0, we find an excellent agreement as shown in 
Fig[TJa). The formation of the spiraling elliptic soliton is due to the effective anisotropic diffrac- 
tion resulting from the OAM, which will be illustrated in the fourth part of the paper, then the 
ellipticity p of the elliptic beam should increase when the critical OAM M c increases, as shown 
in FigOa). In addition, the OAM can strengthen effectively diffraction against self-focusing 
lfl3l . so the critical power P c should increase together with M c when p increases, as shown in 
FigHb). Besides, P c andM t increase when the degree of nonlocality increases, which can also 
be observed in Fig. [2] 




Fig. 1. (color online) Comparison of the beam width of the analytical solution (black 
curves) with that of the numerical simulation (red curves for w x and blue curves for w y ) 
for 8 = 7.5 (a), 5=4 (b), and 8 = 2 (c). 



It can be shown that the solution Eq.© in the strong nonlocality is equivalent to the Gaussian 
complex-variable-function(CVF)-Gaussian soliton, a special case of the C VF-Gaussian solitons 
suggested recently fTT~8l when an arbitrary analytical function / takes the Gaussian function. For 
the limit of the strongly nonlocal nonlinearity, Eq.© and Eq.© can be reduced as <p'| = 

-(£ 2 + C 2 )/(2£ 2 C 2 ),P c |„ = 7TV^ 

use the variable substitutions b 2 = k 2 w 2 /(k 2 + l),c 2 = fC 2 w 2 /(K: 2 — 1), the spiraling elliptic 
soliton (01 can be deduced as i//(x,y,z) = yJP/ (7lbc)f(v) exp [— r 2 / (2w 2 ) — ij3] , where/(u) = 
exp(— V 2 /2) and v = (x + /y)/(K>f)exp(— z'/J). The expression does be the Gaussian CVF- 
Gaussian soliton |[T8l . and the parameter K here is the distribution factor b in Ref. ifTHl . 

3. Analytic stability analysis of the solution 

From Eq.©, we know that P c and a c can be determined when b and c are given. It is also true in 
reverse. When P c and a c are given first, then b and c can be obtained, which are corresponding to 
the stationary point of the potential function Yl(b,c). We use (b s ,c s ) to represent this stationary 
point here. Hence we can study the stability of our analytical soliton solution by determining 
whether (b s ,c s ) is the minimum point of Yl(b, c) or not. For this purpose, we expand n(fr,c) in 
Taylor's series about the stationary point (b s ,c s ) to the second order |fl9l 

(7) 

where jUi = d 2 U/db 2 \ P=Pc ,a=a c , M2 = d 2 U/(dbdc) \ P=Pc>a=C!c , )i 3 = d 2 U/dc 2 \ P=Pc , a=ac , 
Ab = b — b s and Ac = c — c s . We can derive that jii,^ > and \i\ < H1H3 when w m 7^ 0. 
So the stationary point (b s ,c s ) is really the minimum point of the potential function because 
Tl(b, c) — n(b s ,c s ) > in this case. However, we obtain that ji 2 = jiiji^ when w m = 0, and 
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Eq.dZl) is deducedasn(£,c)-n(b. v ,c. v ) « 1/2 [^i {Ab + pL 2 Ac / ^) 2 ]. Thus n(Z>,c) -Il(^,c,) « 
along the particular direction that Ab = — /iiAc /jtti , and next higher order term need to be con- 
sidered in order to judge whether (b s ,c s ) is minimum point. But we can deal with the problem 
in a simpler way. We directly compare the value of the potential function at the stationary point 
(b s ,c s ) and that at the point with a displacement (Ab,Ac) from the stationary point along the 
particular direction that Ab = —\iiAcj[l\, and find that H(b s + Ab,c s + Ac) —H(b s ,c s ) = 0. 
Therefore, for the case of w m = the stationary point (b s ,c s ) is not the minimum point of the 
potential function yet. As a result, we can draw the conclusion that the soliton solutions are 
stable for any degree of nonlocality except for the local case. 

It is known that the spatial profile of the optical beams in the local cubic nonlinear media 
will evolve to a specific circularly symmetric shape, known as the Townes profile |20| that is 
very different from Gaussian profile. So when the degree of nonlocality is weak enough (8 is 
small enough), the Gaussian trial function is not suitable any longer. The deviation between the 
exact solution of Eq.(f2]i and the trial solution Eq.© can bring about the disagreement of the 
beam width obtained from the variational approach with that from the numerical simulation, as 
shown in FigQ2c) when 8 = 2. In fact, the deviation between the variational solution and the 
numerical simulation was discussed in Ref. (T3| for the local case. But, even so, we can find 
that when 5=4, which is out of the region of strong nonlocality, our variational results still 
have a good agreement with the numerical simulations, as shown in FigJTJb). 




Fig. 2. (color online) Critical OAM (a) and critical power (b) as functions of the ellipticity 
p for different degree of nonlocality. 

4. Physical explanation of the formation of spiraling elliptic solitons 

To better understand the formation of such spiraling elliptic solitons, we turn to the analysis of 
the wave vector of the electric field that can be obtained by k = V0 f „/, where 0, o / is the total 
phase of the electric field expressed as fo /(| ,TJ,£) = 0(| , rj, £) +kC, in the physical coordinate 
system. Taking the paraxial beam into consideration, we need only take care of the wave vector 
k around some point on the propagation axis (0, 0, £o)- We therefore can expand {£, , TJ , Q with 
respect to (|,TJ,£) in Taylor's series about (0,0, £o) to the second order, and obtain 

k = {y^ + y^n)^ + {y^n + y^)^+k^, (8) 

where y,- = dj^ =0j r,=ox=^Yji = dj/0l$=o,ii=o,C=&) t/V = 4>T?,0- In me equation above, 
we neglect the terms y^, y^, y^ and y^ because of the fact 11211 that d(j)/d£ <C 
d^/dS, (or d^/drj) for paraxial beams, and take y^ , y 7J = (if y^ , y^ ^ 0, the wave vector 
would have an inclination angle with respect to the £-axis). 

Equation ([8]l tells the fact that the pointing of the vector k at the position , 7] , £) depends 
upon the sign of y^ , y^ , and y^ . For simpleness, we take the projection of k on the , 0, Q- 
plane, representing by k p , into consideration, and the situation for the projection of k on the 



(0, r], Q-plane can be dealt with in the same way. First at the position (<i;,0,£), = 7^§. If 
7^ < 0, we can reach a conclusion that k^ < in the upper half plane (£, > 0) and k^ > in the 
lower half plane (B, < 0). Then k p points downward in the upper half plane and upward in the 
lower half plane, as shown in Fig|2a). As a result, the optical beam will be contracted along the 
4 direction when < 0. When y^% > 0, on the contrary, k p points upward in the upper half 
plane and downward in the lower half plane, as shown in FigfJJb), and the optical beam will be 
expanded along the % direction. For the situation of the position (§,17(7^ 0),£), although the 
presence of the cross-talking term yp„ makes it somewhat complicated, the pointing of k p can 
be determined in the similar way. 




Fig. 3. Schematic of the vector k p at the position (§,0, £) for yes < (a) and 7^ > (b). 
The points we consider are in fact very close to the point (0, 0, £0), but we have magnify 
the extent just for the sake of greater clarity. 



The discussion above about the physical mechanism is based on the analogy of the optical 
beam and the optical pulse. Both of them can be dealt with in Fourier analysis — the tempo- 
ral frequency of the pulse is the analogue of the spatial spectrum of the beam. Therefore, by 
analogy with the phenomenon of the chirp, the time dependence of the phase for the optical 
pulse |221 . the transverse-space dependence of the phase for the optical beam can be referred 
to as "spatial chirp", and the first four terms in Eq. (O are linear spatial chirp terms (the similar 
concept has been introduced in Ref. l23l ). As a result, the physical mechanism for the broad- 
ening (shortening) of the optical pulse ll22l and the expanding (contracting) of the optical beam 
can be understood in this uniform sense. 

On that basis, we discuss the part of the phase due to the OAM expressed as Qoam = 
M sin 2k p ( - % 2 + 77 2 ) / [P (b 2 - c 2 ) ] + 2t, 7]M cos 2/t/3 / [P (b 2 - c 2 ) ] in the phy sical coordinate sys- 
tem, which is corresponding to <Poam = ®XY in the dimensionless coordinate system. Then we 
can obtain the part of the wave vector caused by the OAM 

(oam) 2Msm2kp 2Mcos2kf5 {0 am) 2Msm2kp 2Mcos2/tj3 

*« = ~p(b 2 -c 2 )^ + p{b 2 -c 2 ) 71 ! " = p{b 2 -c 2 ) 7]+ p{b 2 -c 2 y' (9) 

where k^ AM ^ = d^§oAM,k^ AM ^ = d^^oAM- From Eq.© we can find that the contributions of 
OAM to the wavevector are different (asymmetric) in ^-direction and 77 -direction, because the 
signs of the two first terms in Eq. (O are opposite. In other words, OAM can result in an effective 
anisotropic diffraction. It is the effective anisotropic diffraction that leads to the formation of 
elliptic solitons. 

5. Conclusion 

We have obtained spiraling elliptic solitons in nonlocal nonlinear media without anisotropy by 
use of the variational approach. The formation of such solitons is due to an effective anisotropic 
diffraction resulting from the orbital angular momentum. We show that this class of solitons are 



stable for any degree of nonlocality except the local case. Our approximate analytical results 
have been confirmed by direct numerical simulations of the NNLSE. 

Acknowledgments 

This research was supported by the National Natural Science Foundation of China (Grant Nos. 
11074080 and 10904041), the Specialized Research Fund for the Doctoral Program of Higher 
Education (Grant No. 20094407110008), and the Natural Science Foundation of Guangdong 
Province of China (Grant No. 1 1 5 1 063 1 1 00001 7). 



